Magnetic interactions in transition-metal doped ZnO: An ab-initio study 
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We calculate the nature of magnetic interactions in transition-metal doped ZnO using the local 
spin density approximation and LSDA+ U method of density functional theory. We investigate the 
following four cases: (i) single transition metal ion types (Cr, Mn, Fe, Co, Ni and Cu) substituted 
at Zn sites, (ii) substitutional magnetic transition metal ions combined with additional Cu and Li 
dopants, (iii) substitutional magnetic transition metal ions combined with oxygen vacancies and (iv) 
pairs of magnetic ion types (Co and Fe, Co and Mn, etc.). Extensive convergence tests indicate that 
the calculated magnetic ground state is unusually sensitive to the k-point mesh and energy cut-off, 
the details of the geometry optimizations and the choice of the exchange-correlation functional. 
We find that ferromagnetic coupling is sometimes favorable for single type substitutional transition 
metal ions within the local spin density approximation. However, the nature of magnetic interactions 
changes when correlations on the transition-metal ion are treated within the more realistic LSDA + 
U method, often disfavoring the ferromagnetic state. The magnetic configuration is sensitive to the 
detailed arrangement of the ions and the amount of lattice relaxation, except in the case of oxygen 
vacancies when an antiferromagnetic state is always favored. 
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I. INTRODUCTION 

Dilute magnetic semiconductors (DMSs), obtained by 
partial replacement of the cations in conventional semi- 
conductors by magnetic transition-metal ions, are of cur- 
rent interest as potential semiconductor-compatible mag- 
netic components for spintronic applications 1 . Early 
studies of DMSs focussed on II- VI semiconductor hosts 
such as CdSe and ZnTe 2 ; since many transition met- 
als adopt divalent ionic states, they therefore substitute 
readily for divalent cations such as Zn 2+ or Cd 2+ . How- 
ever, most II- VI DMSs are antiferromagnetic or have 
very low ferromagnetic Curie temperatures 3 -, rendering 
them unattractive for applications. More recently, robust 
ferromagnetism was observed in DMSs based on III-V 
semiconductors such as Mn-doped GaAs, in which Curie 
temperatures of ~150 K have been achieved^. Here, 
the ferromagnetic coupling between the localized Mn 
magnetic moments is believed to be mediated by itiner- 
ant carriers (holes) which are introduced when divalent 
Mn ions replace trivalent galliums 5 ^ (so-called carrier- 
mediated ferromagnetism) . Calculations for other mate- 
rials, based on such a carrier-mediated mechanism with 
a high concentration (~5 %) of holes, predicted above- 
room-temperature ferromagnetism for Mn-doped ZnO 
and GaN 5 . 

In addition to the technological appeal of room- 
temperature ferromagnetism in ZnO-based DMSs, ZnO 
offers other desirable features as a semiconductor host. 
It has a direct wide band gap of 3.3 eV— ,and therefore 
finds widespread use in the optoelectronics industry 8,9,10 . 
Its strong piezoelectricity 11 is exploited in a variety of 
transducer applications^ and have possible application 



in polarization field effect transistors 13 :. And long spin 
coherence times, with potential spintronic applications, 
have recently been reported at room temperature in n- 
type ZnQii Thus, if it could be achieved, ferromagnetic 
ZnO would be a highly multifunctional material with co- 
existing (and possibly coupled) magnetic, piezoelectric, 
optical and semiconducting properties. 

The predictions of high-temperature ferromagnetism- 5 
spawned a large number of experimen- 
tal 15 : 16 . 17 ! 18 . 19 : 20 . 21 ! 22 ! 23 ! 24 ! 25 and computa- 

tional 26 ! 27 ! 28 ! 29 ' 30 ! 31 ' 32 studies of transition-metal 
(TM)-doped ZnO. The reported experimental values 
of Curie temperature and magnetization show a large 
distribution, suggesting that the system is sensitive 
to preparation methods, measurement techniques, 
substrate choice, etc. For example, in the case of 
Mn-doped ZnO, while some experiments report above 
room temperature ferromagnetismSSi 3 ^, others report 
a low ferromagnetic ordering temperatureiS^i^ or a 
spin-glass or paramagnetic behavioaiSiS 2 * 3 ^. Likewise, 
in Co-doped ZnO, there are reports of giant magnetic 
moments (~ 6.1 /is/Co)2i, high ferromagnetic ordering 
temperatures with moments of ~ 1-3 ^^Agj^LSiiS and 
in some cases no ferromagnetic behavior at gj^^^J^^ 
Interestingly, the reported computational results show 
a similar spread of values; this arises in part from 
different physical approximations (choice of exchange- 
correlation functional, inclusion or omission of structural 
relaxations), but also from the unusual sensitivity of 
the magnetic interactions to the convergence quality 
of the computations. (For a detailed review see Ref. 
l24j) . As a result, in spite of the flurry of experimental 
and theoretical work, no definite conclusions have been 



reached regarding the nature and origin of the magnetic 
interactions in this system. 

In this work, we report our results of a systematic com- 
putational study of transition- metal-doped ZnO, with the 
goal of understanding the nature and origin of the mag- 
netic interactions. Our emphasis is on extracting trends 
along the 3d transition metal series, and so we first ex- 
plore the effects of substitutional doping with a range of 
transition metals (Cr, Mn, Fe, Co, Ni and Cu). Next, we 
calculate the effects of additional dopant or vacancy im- 
purities, and finally we calculate the preferred magnetic 
configurations of pairs of different magnetic ion types. 
Our main finding is that the ferromagnetic state is gen- 
erally disfavored in the most realistic calculations, and 
it is not strongly stabilized by common defects likely 
to be found in as-grown ZnO, or easily incorporated as 
dopants. 

The remainder of this paper is organized as follows. 
In the next section, we describe our computational and 
system details and outline the unusually demanding con- 
vergence behavior that we find for TM-doped ZnO. In 
Sec. 11111 we present our calculated trends across the 3d 
series for single-type transition metal doping within the 
local spin density approximation (LSDA) and the LSDA 
+U method. In Sec. II VI we investigate the influence of 
defects across the 3d series by (i) co-doping with Cu or 
Li and (ii) incorporation of oxygen additional oxygen va- 
cancies. In Sec. |V] we present our results for double 
doping (simultaneous doping of two different TM ions) 
in (Zn,(Co,Fe))0 and (Zn,(Mn,Co))0. Finally, we sum- 
marize our results in Sec. I VII 



II. TECHNICALITIES 



A. System details 



In order to achieve realistic experimental dopant con- 
centrations (~ 10 - 30%), we used a periodic 2x2x2 
wurtzite supercell of ZnO which consists of 32 atoms in 
a unit cell. Substitution of two Zn atoms by transition 
metal ions then gives a dopant concentration of 12.5% 
and allows for calculation of the relative energies of ferro- 
magnetic (FM) and antiferromagnetic (AFM) orderings. 
We explored two spatial arrangements, near (transition 
metal atoms separated by one oxygen atom) and far (TM 
atoms separated by -O-Zn-O- ), as shown Figure [21 In 
each case, we used the energy difference between FM 
and AFM ordering, AE = Eafm-^fm, as an indica- 
tor of the magnetic stability (a positive AE implies that 
FM is favorable). Equating AE with the thermal energy 
fcflT suggests room temperature ferromagnetism should 
be achieved for AE values larger than ~ 30 meV. 
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FIG. 1: Plot of Eafm-Efm for Zno.875Coo.125O as a function 
of k-point grid, and for three different plane- wave energy cut- 
offs (ecut). 



B. Method details 

Our total energy and electronic structure calculations 
were performed using the projector augmented wave 
(PAW) 41 formalism of density functional theory as im- 
plemented in the VASP packago 42 i 43 i 44 . We used the de- 
fault VASP PAW potentials with 12 valence electrons for 
Zn (3d 10 4s 2 ), 6 for O {2p 4 2s 2 ) and 2+n (3d n 4s 2 ; n = 3 
to 9) for the transition metals. We used a well converged 
energy cut-off of 550 eV for plane-wave expansion of the 
PAW's, a 4 x 4 x 3 Gamma centered k-point grid, and 
the tetrahedron method with Blochl 4 ^ corrections for the 
Brillouin Zone integrations. Note that this energy cut-off 
and k-point sampling are unusually high; our convergence 
tests (Figure ^ indicate that qualitatively incorrect mag- 
netic behavior often occurs for lower values. Similarly 
rigorous convergence parameters have been shown to be 
required for accurate calculations for Co-doped Ti02^- 
For all geometry optimizations, we kept the volume of 
the supercell fixed to the experimental volume (a=6.50A 
and c=10.41 A£) and relaxed all the internal coordinates 
until the Hellman-Feynman forces were less than 10~ 3 
eV/ A. 

We approximated the exchange-correlation functional 
with both the local spin density approximation (LSDA) 
and the fully localized limit of the LSDA+ U method 46 . 
Although widely used and well established for many 
properties, the LSDA is well-known to yield incorrect 
behavior for strongly correlated magnetic systems, often 
predicting half-metallic, low spin states for systems which 
are actually insulating and high spin. For example, most 
TM monoxides are wide band gap antiferromagnetic in- 
sulatora 47 ! 48 ! 49 ! 50 ! 51 , however the LSDA finds them to be 
either FM metals (FeO and CoO) or small-gap semi- 
conductors (MnO and NiO)£2. The LSDA+[/ method 
extends the LSDA by explicitly adding the on-site d-d 
Coulomb interaction, U, and the on-site exchange inter- 
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action, J, to the LSDA Hamiltonian, and usually gives 
improved results for magnetic insulators. Here we use 
typical values of U = 4.5 eV and J — 0.5 eV on all tran- 
sition metals^*^^; although a small variation in U and 
J across the series is expected, our choice of constant 
values permits a more straightforward comparison. We 
do not modify the LSDA ZnO electronic structure, which 
therefore shows the usual LSDA underestimation of the 
band gap, and underestimation of the Zn d state energies. 




(a) (b) 

FIG. 2: (Color online) Wurtzite supercell of ZnO with two TM 
ions in the (a) near and (b) far configuration. The (yellow) 
light-shaded spheres are the TM ions, the large (cyan) dark- 
shaded spheres the Zn ions, and the small (red) black spheres 
the O ions. 



III. MAGNETIC INTERACTIONS IN DOPED 
ZNO 

A. LSDA results 

We begin by calculating the total energies of 12.5% 
TM-substituted ZnO within the local spin density ap- 
proximation (LSDA). We reiterate that we do not expect 
the LSDA to give accurate magnetic behavior for this 
system, but present the results as a baseline for com- 
parison with our LSDA+U results in the next section. 
Keeping the volume fixed at the ZnO experimental vol- 
ume, we first relax all internal coordinates for near and 
far arrangements of the TM ions, and for both AFM and 
FM orderings. In figure[3]we show our calculated energy 
differences between the near and far configurations; it is 
clear that in all cases (except Mn, for which the energy 
difference is negligible) it is favorable for the TM ions to 
cluster together. 

Figures EI a ) and (b) show our calculated magnetic 
energy differences (Eafm-^fm) for our range of TM 
dopants in the near and far spatial arrangements. First 
we note that our calculated magnetic orderings without 
ionic relaxations, which vary considerably both with the 
TM type and with the spatial arrangement, are consis- 
tent with earlier calculations in the literaturei£*2L£&. We 
see, however, that the strength and sign of the mag- 
netic interaction is highly sensitive to the ionic relax- 
ation, which changes both the distance between the TM 
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FIG. 3: LSDA energy differences (AE = E near - E far ) for 
substitutional TM ions. The negative energy difference indi- 
cates that the TM ions prefer to be in a near spatial config- 
uration. 

ions (by up to ~ 0.25 A) and the TM-O-TM angle (by 
up to ± 5%). We find no clear trends in the magni- 
tude of relaxation, nor any direct correlation between the 
relaxation and the change in the magnetic interactions, 
across the TM series. We attribute the sensitivity in mag- 
netic ordering to a subtle competition between antiferro- 
magnetic superexchange (favored by 180° bond angles), 
ferromagnetic superexchange (favored by 90° TM-O-TM 
bond angles) and AFM direct exchange (favored by short 
TM-TM distances) 57 i 58 i 59 . Our results indicate that ionic 
relaxations must always be included in calculations to ob- 
tain meaningful results, and that even the generalization 
of relaxed positions from one TM ion to another should 
be applied with cautioniS^. For the most realistic con- 
figurations {near, with relaxations), the LSDA suggests 
that doping with Cr, Fe, Co, Ni and Cu should lead to a 
ferromagnetic ground state. 



Next (Figure|3J) we compare our calculated total densi- 
ties of states (DOS) and TM 3c? projected local densities 
of states (PDOS) for our range of TM ions in ZnO. In all 
cases the DOS for the FM arrangement is shown. The 
total DOS is represented by the gray shaded region while 
the black shaded regions represent the 3d states of the 
TM impurities. The majority (f) spin states are plot- 
ted along the negative x direction, and the minority (J,) 
states are plotted along the positive x direction. The en- 
ergies are reported relative to the Fermi energy (E/=0). 
For comparison, the DOS of undoped ZnO is also shown, 
with the Fermi energy set to the top of the valence band. 

Comparison with the ZnO DOS identifies the broad 
(~ 4-5 eV) band just below the Fermi energy as derived 
largely from the O 2p states, with the narrow Zn 3d band 
just below and slightly overlapping with the O bands. 
The bottom of the conduction band is composed largely 
of Zn 4s states, and the band gap (^0.78 eV) shows the 
usual LSDA underestimation. 

In all cases, the exchange-split TM 3d states form fairly 
narrow bands in the region of the gap. The exchange 
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FIG. 4: LSDA energy differences (AE = Eafm — Efm) for substitutional TM ions in ZnO for (a) near and (b) far spatial 
configurations. In both (a) and (b), the filled bars represent the case when the ionic coordinates are not relaxed and the unfilled 
bars represent the case when the ionic coordinates are relaxed. Lines at AE = indicate that FM and AFM orderings are 
equivalent in energy. 



splitting (~2 eV) is consistently larger than the crystal 
field splitting (~0.5 eV) which splits each spin manifold 
into lower energy doubly degenerate e and higher en- 
ergy triply degenerate t states. We observe the following 
trends, which are consistent with the crystal chemistry 
of 3d transition metal oxides: 

• The energy of the TM 3d states relative to the top 
of the valence band shifts down in energy on moving 
right across the 3d series (from Cr to Cu). As a con- 
sequence, the Cr states are somewhat hybridized 
with the bottom of the conduction band (this is 
likely a result of the LSDA underestimation of the 
ZnO band gap), the Fe and Co states are largely 
mid-gap, and the Cu states are hybridized with the 
top of the valence band. 

• The calculated DOSs (Figure[5J are consistent with 
the Hund's rule predictions for isolated TM ions, 
with only majority-spin (f ) states occupied for Cr 
and Mn, and occupation of the minority-spin (J.) 
states increasing Fe to Co to Ni to Cu. Also the 
exchange splitting decreases on moving to the right 
across the 3c? series. 

• The (d 5 )T Mn 2 + ion, and (d 5 )T (d 2 )^ Co 2 + ions lead 
to insulating behavior; the other ions have partially 
filled d manifolds and a finite density of states at 
the Fermi level. 

• The variation in magnetic moments across the 3d 
series reflects the decrease in exchange splitting, 



and the increase in number of electrons, with the 
local magnetic moments (obtained by integrating 
up to a Wigner- Switz radius of 1.4 A)S& on the TM 
ion as follows (in units of /is): 3.04 (Cr), 4.31 (Mn), 
3.51 (Fe), 2.44 (Co), 1.49 (Ni) and 0.55 (Cu). 



These LSDA results are consistent with earlier LSDA 
calculationsiZi2L2£, and also with some experimental re- 
ports^. Since the LSDA is known to underestimate the 
band gaps and exchange splittings in magnetic systems, 
we withold a detailed analysis of the band structures until 
the next section, where we repeat our calculations with 
the inclusion of correlations at the level of the LSDA + 
U method. 



B. LSDA + U 

Next we repeat the suite of calculations described 
above, with the correlation on the 3d TM ions treated 
within the LSDA +U scheme 46 . Here we present results 
obtained using typical values of U = 4.5 eV and J = 0.5 
eV; in the appendix we discuss the U -dependence of our 
calculated properties. 



First, in Figure H3 we present our calculated LSDA+U 
DOSs of TM-doped ZnO. Although all the trends across 
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FIG. 5: DOS and PDOS of TM-doped ZnO with the TM atoms in the near spatial configuration calculated within the LSDA. 
The black shaded regions show the TM d states, and the gray shaded regions show the total DOS. For clarity, the TM d states 
are scaled by a factor of two. Also shown (right panel) is the DOS for undoped ZnO. 
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FIG. 6: DOS and PDOS of TM-doped ZnO with TM atoms in the near spatial configuration calculated within the LSDA + U 
method. The black shaded regions show the TM d states, and the gray shaded regions show the total DOS. For clarity, the 
TM d states are scaled by a factor of two. Also shown (right panel) is the DOS for undoped ZnO. 



the series remain the same as in the LSDA, the addition 
of the Hubbard U term, as expected, causes an increased 
splitting between the filled and empty orbitals in the TM 
d manifold. The exchange splitting is now increased to 
~4 eV, and the crystal- field splitting increases to ~1 eV. 
The TM d states, which were localized in the gap in the 
LSDA, now shift down in energy and hybridize strongly 
with the O 2p states. Comparing the DOS of Fe and Ni, 
we observe a transition from half-metallic to insulating 
when going from the LSDA to LSDA +TJ. This is directly 
reflected in the magnetic interation strength. For e.g., in 
Fe and Ni, the long-range interaction vanishes and the 
interaction strength decreases. We also calculated the 
magnetic moments on the TM ion as follows (in units 



of fj, B ): 3.40 (Cr), 4.41 (Mn), 3.54 (Fe), 2.56 (Co), 1.63 
(Ni) and 0.65 (Cu). The magnetic moments are slightly 
larger compared to the LSDA values. Next we compare 
the total energies of the FM and AFM magnetic orderings 
in the near and far spatial configurations. In all cases, 
all the ionic coordinates are relaxed within the LSDA +U 
scheme. Figures Ufa) and (b) show the energy difference 
between the FM and AFM states, (Eafm — Efm), for 
near and far arrangements for the various TM ions. As 
in the LSDA case, AE shows a strong dependence on 
the spatial arrangement of the dopant ions, as well as on 
the extent of ionic relaxations included in the calculation. 
In all cases, the sign and magnitude of the interactions 
are strikingly different from those obtained within the 
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FIG. 7: LSDA +U energy difference, AE=E afm — Efm, for substitutional TM ions in ZnO represented for (a) near and (b) 
far spatial configurations. The filled bars show the energies when the ions are not relaxed and the shaded bars give the values 
for when relaxations are included. Lines at AE=0 indicate that FM and AFM orderings are equivalent in energy. 



LSDA, with a general trend to reduced ferromagnetism. 
For the energetically favorable near case, only Ni and Cu 
dopants show a tendency to ferromagnetic ordering when 
ionic relaxations are included in the calculation. 



IV. INFLUENCE OF DEFECTS 



Possible p-type dopants; Cu and Li with TM in 
ZnO 



Given that our calculations containing single types of 
subsitutional dopants are unable to reproduce the exper- 
imentally reported ferromagnetism, we next search for 
other impurities that could mediate the ferromagnetic 
ordering. Since the original predictions of high Curie 
temperature in Mn-doped ZnO^ assumed high hole con- 
centrations, we first include additional dopants which are 
likely to introduce p-type carriers. Some experimental 61 
and theoretical studies^ 6 - have suggested that Cu could 
provide holes when doped into ZnO. Cu has an electronic 
configuration of 3d 9 4s 2 , and has a tendency to form Cu 1+ 
ions with a 3d 10 4s° electron configuration, and an ionic 
radius (0.60 A) close to that of Zn 2+ . Thus it is rea- 
sonable to assume that substitution of Zn with Cu could 
yield a Cu 1+ configuration with an associated hole. In- 
deed, earlier LSDA calculations on Cu-doped ZnO found 
a band structure consistent with this model, and an en- 
hanced tendency to ferromagnetism^. Lithium has also 
been suggested as a possible p-type dopant in ZnQ£2. Li 
has an electronic configuration of Is 2 , 2s 1 , and, like Cu, 



tends to form a singly charged cation (ls 2 2s°) with an 
ionic radius (0.59 A) close to that of Zn 2+ . In this section, 
we calculate the effects of co-doping TM-doped ZnO with 
Cu or Li, again within the LSDA+U method, using the 
same convergence and U and J parameters as described 
above. Our calculated AE values are shown in figure |S1 
In all cases the concentration of the Cu or Li was 6.25% 
(one additional dopant in 16 cations), and the additional 
dopant was placed far away from the TM ion in ZnO. It 
is clear that, in contrast with earlier LSDA calculations, 
co-doping with Cu within the LSDA+U method does not 
strongly enhance the tendency towards ferromagnetism, 
yielding an antiferromagnetic ground state in most cases, 
and a weakly ferromagnetic state for Cr and Ni. Co- 
doping with Li, however, is more favorable, and in some 
cases should lead to above-room-temperature ferromag- 
netism. This is particularly intriguing in light of recent 
reports of ferroelectricity in Li-doped Zn Q 63 i 64 , suggest- 
ing the possibility of engineering a multiferroic material 
with simultaneous ferromagnetism and ferroelectricity. 
However, none of our calculated Li-doped materials is 
insulating, a requirement for ferroelectricity. 



B. O vacancies: Zn(TM,V )0 

Next, we combine substitutional TM ions with oxy- 
gen vacancies, which are believed to be the most com- 
mon native defects in ZnO 62-65,66 . There is some ex- 
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FIG. 8: Eafm-Efm for TM-doped ZnO with 6.25% Cu/Li. 



pcrimental evidence associating the presence of oxygen 
vacancies with the existence of ferromagnetism in TM- 
doped ZnO. For example, Venkatesani& et al. reported a 
correlation between the magnitude of Co magnetic mo- 
ments and the oxygen partial pressure during annealing, 
with higher oxygen partial pressure reducing the amount 
of magnetization. In addition, Kittilsved et a/ 23 report 
the observation of room-temperature ferromagnetism in 
Co-doped ZnO nanocrystals with oxygen-containing sur- 
faces, and propose that oxygen vacancies mediate the 
magnetic interactions. A couple of experiments also ob- 
serve a change in the magnetization in samples annealed 
in reduced atmospheric pressures and they attribute it 
to oxygen vacancies^. As a result of this reported cor- 
respondence between the presence of oxygen vacancies 
and the existence of ferromagnetism, a model has been 
proposed in which ferromagnetism is mediated by carri- 
ers in a spin-split impurity band derived from extended 
donor orbitalsiS^. The validity of the model rests on 
the formation of an oxygen-vacancy-derived donor im- 
purity band close to the ZnO conduction band edge, 
which hybridizes with the spin-polarized TM 3c? band. 
If the Fermi energy lies within this hybrid oxygen va- 
cancy/TM 3d spin-polarized band, a carrier-mediated 
ferromagnetism should be favorable. Here we calculate 
the relative energies of the oxygen impurity levels and the 
transition metal 3d states across the TM series, in order 
to investigate the applicability of this model. Indeed, by 
simple charge neutrality arugments, removal of a neutral 
oxygen atom should leave two unbonded electrons associ- 
ated with the vacancy, provided that the Zn remains in a 
Zn 2+ state. We note, however, that recent first-principles 
calculations- 8 found oxygen vacancies to be deep donors 
in ZnO, and therefore unavailable for carrier mediation. 

We introduce a single O vacancy as far as possible from 
the TM ions which are again in both the near and far 
configurations. The energy difference {Eafm — -^fm) is 
plotted in Figure [5] We see that the AFM state is stable 
in all cases, consistent with previous computations for 
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FIG. 9: Eafm — Efm for the range of TM ions in the pres- 
ence of an oxygen vacancy. The shaded and unshaded bars 
represent the cases when the TM ions are near and far from 
each other respectively. In both cases the oxygen vacancy is 
placed as far as possible from the TM ions. 



Co-doped ZnO with a neutral oxygen vacanc; 



,17,29,30 



In Figure HTJI we show our calculated DOSs for the en- 
tire series. Note that the defect concentration used in 
our calculations is higher than likely attainable experi- 
mentally (~l-2%), but that the relative band alignments 
will be largely unchanged by the increased concentration. 
We see a shift of the Fermi energies relative to those in 
the absence of oxygen vacancies, but no other striking 
changes in band structure. In particular, it is clear that 
the Fermi level does not lie in a hybridized TM 3d-0 2p 
impurity band for any case. Therefore our band struc- 
tures are not consistent with those required to mediate 
ferromagnetism within the model proposed above; they 
are however consistent with our calculated AFM ground 
states. 

The calculated magnetic moments (in units of /zb) for 
each TM ion in the presence of the oxygen vacancy are 
3.54 (Cr), 4.36 (Mn), 3.66 (Fe), 2.64 (Co) and 1.68 (Ni), 
largely unchanged from the LSDA+U values without the 
oxygen vacancy. This indicates that the two electrons 
from the oxygen vacancy are localised and do not influ- 
ence the occupancy of the TM d states. Finally, we point 
out that we have considered only neutral oxygen vacan- 
cies here. A recent first principles study suggests that 
positively charged oxygen vacancies might be more suc- 
cessful for mediating ferromagnetism in Co-doped ZnO 30 . 



C. Other defects 

Next, we study the effect of a range of experimentally 
plausible defects - Zn vacancies (Vzn), octahedral Zn 
interstitials (Zn,), TM interstitials (Co,;) and Li intersti- 
tials (Lii) - on the magnetic interactions. Rather than 
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FIG. 10: DOS and PDOS for TM-doped ZnO (near spatial arrangement and FM state) in the presence of oxygen vacancies. 
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FIG. 11: Eafm — Efm for ZnO doped with Co and a range 
of other defects. 



scanning the entire range of transition metal dopants, 
we use Co-doped ZnO as our test system, since it is the 
most widely studied both experimentallyi£ii£i22iSL3i and 
theoreticallyiL22iiffii££i£i. Figure ITT1 shows our calculated 
AE = Eafm — Efm for the near spatial arrangement, 
with the vacancies placed as far as possible from the Co 
ions. Our previously reported values for substitutional 
Li (Lizn) and oxygen vacancies (Vo) are also shown for 
completeness. 



Perhaps the most important result is that, as in the case 
of substituional Li reported above, interstitial Li also sta- 



bilizes the ferromagnetic state. This is significant since 
Li impurities are likely to be incorporated both in substi- 
tutional and interstitial sites during growth, and avoid- 
ing competing antiferromagnetic interactions is desirable. 
We also see that Zn vacancies are favorable for ferro- 
magnetism, consistent with previous reports in the liter- 
atur e) 1 ^ 29 ' 30 . The other defects studied do not stabilize 
the ferromagnetic state. 



V. SIMULTANEOUS DOPING OF CO WITH 
OTHER TM IONS 

Finally, we calculate the properties of ZnO doped 
with both Co and an additional TM ion (Mn or Fe). 
There has been some experimental work on such co- 
doped ZnO-based DMSs. Cho£& et al. reported room- 
temperature magnetism with a saturation magnetization 
of 5.4 emu/g (corresponding to ~ lfis for CoFe pair) for 
Zni_ x (Feo.5,Coo.5)a;0 (x=0.15) films fabricated by reac- 
tive magnetron co-sputtering. Han- 1 et al. also reported 
room temperature magnetism for Zni_ a; (Fei_j ) ,Cuy) a; 
bulk samples with To ^550K. On the theory side, Park 
and Min2& calculated the electronic structures and mag- 
netic properties of (Fe,Co) and (Fe,Cu) doped ZnO at 
12.5% doping. They reported a tendency to form FM Fe- 
O-Cu clusters, and argued that ferromagnetism should 
arise from double exchange, whereas the absence of clus- 
tering found for (Fe,Co)-doped ZnO would require a dif- 
ferent mechanism for ferromagnetism. 

Here we explore the possibility of ferrimagnetic (FiM) 
ordering, in which the magnetic moments of one TM 
ion type are antiparallel to those of the other TM ion 
type, but a net magnetization arises from incomplete 



9 





FiM 



FM 



AFM 



FIG. 12: (Color online) A schematic of the different magnetic 
orderings investigated in this work. The black arrows rep- 
resent the magnetic moments of the Co ions while the blue 
(gray) arrow represent the magnetic moment of the other TM 
ion (Mn or Fe). 




> 

>> 

SIJ 
u 

s 




FIG. 13: DOSs and PDOSs of 3d states for ferrimagnetic 
Zno.875(Coo.5Mno. 5)0.1250 and Zno.875(Coo.sFeo. 5)0.1250. The 
blue shaded region represents the d states of the Co ions while 
the red shaded region represents the d states of the TM ions. 



cancellation of the magnetic moments. For each sys- 
tem (Zn(Co,Fe)0, Zn(Co,Mn)0) we calculate the total 
energies of the FM, AFM and FiM magnetic orderings 
as shown in Figure El with TM ions of different type 
in the near arrangement. We use a 4x4x2 wurtzite su- 
percell with 64 atoms (double the size of the previous 
studies), which allows the permutations shown in Fig- 
ure El a t a TM ion concentration of 12.5%. Within this 
constrained spatial arrangement, we find that in both 
systems the FiM ordering is the most stable arrange- 
ment by 200 meV for (Zn,(Co,Fe))0, and by 57 meV for 
(Zn,(Co,Mn))0; the corresponding magnetic moments 
per supercell are 1.11 \ib for (Zn,(Co,Fe))0 and 1.98 
\ib for (Zn,(Co,Fe))0 which agrees well with the ex- 
perimentally reported values™. While ferrimagnetism 
overcomes the problem of cancellation of magnetic mo- 
ments by superexchange-driven antiparallel alignment, it 
of course requires that ions of different types cluster to- 
gether while those of the same type remain distant; this 
might be difficult to achieve experimentally. 

For completeness, we show the DOSs of FiM 
Zno. 8 75(Coo. 5 Mo.5)o.i250 (M = Mn, Fe) in Figure [TJ 
The (Zn,(Mn,Co)) DOS is almost an exact superposition 
of the separate (Zn,Mn)0 and (Zn,Co)0 DOSs. Inter- 
estingly, in the (Zn,(Fe,Co))0 system, we get a metallic 
state because of the narrow band spin polarized impu- 
rity band at the Fermi level unlike either (Zn,Fe)0 or 
(Zn,Co)0. The Coey model described in the previous 
section could be appropriate in this picture where the 
donor impurity band of one the transition metal atom 
(in this case Fe) lies close to the conduction band edge 
and overlaps with the Fermi level. Thus we have a stable 
ferrimagnetic state. 



VI. SUMMARY 



In summary, we have performed a systematic study of 
the magnetic behavior of TM-doped ZnO for a range of 
TM ions and defects. Our main result is the absence, in 
general, of a tendency for pairs of TM ions substituted 
for Zn to order ferromagnetically; in most cases AFM 
ordering is more favorable. Ferromagnetic ordering of 
TM ions is not induced by the addition of substitutional 
Cu impurities nor by oxygen vacancies. Incorporation 
of interstitial or subsitutional Li is favorable for ferro- 
magnetism, as are Zn vacancies. On a technical note, 
we find that the calculated magnetic behavior is strongly 
dependent both on the computational details (with ferro- 
magnetism disfavored by improved convergence) and on 
the choice of exchange-correlation functional (with ferro- 
magnetism disfavored by the more appropriate LSDA+U 
method). This observation explains the large spread of 
computational results in the literature. 
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APPENDIX: DEPENDENCE OF PROPERTIES 
ON CHOICE OF U: THE EXAMPLE OF 
CO-DOPED ZNO 

Here we illustrate the sensitivity of our results to the 
magnitude of the U parameter with calculations of the 
energy difference, Eafm ~ Efm and total and orbital- 
resolved densities of states, as a function of U in the 
(Zn,Co)0 system. The results are plotted in Figure PHI 
Figure ITU (upper panel) shows E AFM — E FM for the 
LSDA case (with U=0, (i)), for U on the Co d states 
ranging from 4 to 6 eV (points (ii), (hi) and (iv)) and 
with a U of 4 eV on both the Co and Zn d states (v). 
Two conclusions can be drawn from the plot: First, the 
sign of the magnetic interaction changes from LSDA to 
LSDA+ U in the near case making the AFM ordering 
more stable. The AFM remains stable for a range of 
U values. Second, the U on the Co ion dominates the 
magnetic interaction; the energy differences in cases (ii) 
and (v) are very close, indicating that correlations on 
the Zn d states do not significantly affect the magnetic 
interactions. FigurelbH flower panel) shows the DOS and 
Co d PDOS in Co-doped ZnO for a range of U values. 
We see that, as expected, the exchange splitting between 
occupied majority states, and unoccupied minority states 
increases with U; as a consequence the majority occupied 
states move down into the valence band and hybridize 
more strongly with the O 2p states as U is increased. 
Adding a U of 4.5 eV on the Zn also lowers the energy 
of the Zn d states. 
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FIG. 14: Upper panel: Total energy difference (EafM-fu) 
for 12.5% Co in ZnO for various U values on the Co d states. 
The last bar gives this energy difference with U = 4 eV on 
both Zn and Co 3d states. The gray shaded and the black 
shaded bars show the energy differences for the near and far 
cases. Lower panel: Total DOS of (Zn,Co)0 and PDOS of 
Co d states for a range of U values. The light shaded re- 
gions represent the total DOS while the black shaded regions 
correspond to the 3d states of Co ion. 
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